One-dimensional compression of Bose-Einstein condensates by laser-induced 

dipole-dipole interactions 
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We consider a trapped cigar-shaped atomic Bose-Einstein condensate irradiated by a single far-off 
resonance laser polarized along the cigar axis. The resulting laser induced dipole-dipole interactions 
between the atoms significantly change size of the condensate, and can even cause its self-trapping. 
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The ability to alter the interatomic potential of weakly interacting Bose-Einstein condensates (BECs) by applying 
electromagnetic fields affords a high degree of control over the Hamiltonian and thus the possibility to engineer the 
macroscopic properties of such many-body systems ^ . The experiment of Inouye et al |^ has demonstrated how the 
s-wave scattering length can be changed by magnetic fields via a Feshbach resonance. 

We have recently proposed a different avenue: the use of off-resonant lasers to induce long-range dipole-dipole 
interactions in atomic BECs These interactions can drastically modify a BEG, causing its self-trapping laser 

induced self-" gravity" ||] (for certain laser-beam configurations), "supersolid" structures js), and pecuhar excitation 
spectra At the same time, off-resonant lasers may undergo "superradiant" Rayleigh scattering from the 
condensate, concurrently with collective atomic recoil (GARL) It is important to draw a clear distinction between 
such effects and those of laser-induced dipole-dipole forces. The purpose of this paper is to suggest a simple geometry 
in which the compression of a trapped condensate and its possible self-trapping are unambiguously caused by laser- 
induced dipole-dipole interactions, whereas superradiant Rayleigh scattering or GARL are absent. 
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FIG. 1. The laser beam and condensate geometry. By choosing the polarization to be along the long axis of the condensate 
superradiant effects are suppressed. 

A trapped cigar-shaped BEG that is tightly confined in the radial x-y plane is irradiated by a plane-wave laser at a 
wavelength larger than the radial size of the condensate. The laser polarization is chosen to be along the long z-axis 
of the condensate in order to suppress "superradiant" Rayleigh scattering or GARL Q that are forbidden in the 
direction of the field polarization. The interatomic dipole-dipole potential induced by far-off resonance electromagnetic 
radiation of intensity /, wave- vector q = (/y (along the y-axis), and polarization e = z (along the z-axis) is P| 



C/(r) 



a^{q)Vzziq,r) cos(gy). 



(1) 



Here r is the interatomic axis, a(q) the isotropic, dynamic, polarizability of the atoms at frequency cq, and Vzz is the 
appropriate component of the retarded dipole-dipole interaction tensor 



Vzz = 4 



(l — 3cos^(0)) (cosqr + singr) — sin^ {9)q^r^ cos qr 



(2) 



where 6 is the angle between the interatomic axis and the z-axis. We note that the far-zone {qr 3> 1) behavior of the 
dipole-dipole potential (|l|) along the z-axis direction is proportional to — sin{qr) / (qr)'^ . 

A mean-field description of a BEG with dipole-dipole forces QJl^,|ll| can be accomplished through the Gross- 
Pitaevskii equation [n2[ for the condensate order parameter 



ih 



'dt ~ S^*^ 



-Ht. 



(3) 
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The mean-field energy functional 

^^tot = -f^kin + ^^ho + ^^dd + ^^s (4) 

contains contributions from the one-body energies and the interatomic potentials: (a) the kinetic energy 

i?ki„ = J {h^/2m)\V^^d^r; (5) 

(b) harmonic trapping 

Hho = Jvi,on{r)d\, (6) 
where n(r) is the atomic number density; (c) the dipole-dipole potential 

Hdd = (1/2) / n(r)n(r') C/(r - r') d^r d^r' ; (7) 



(d) the mean-field energy due to the very short-range (r ^) van der Waals interaction, which will be treated, as is 
usual, within the delta function pseudo-potential approximation, 

= (l/2)(47ra?iV"^) j n{r)d^r ; (8) 

where a is the s-wave scattering length. 

The integrations involved in the energy functionals are simpler in momentum space, for which H^d = (l/2)(27r)'^/^ 
/ C/(k) n(k) «(— k) d^k. The Fourier transform of the dipole-dipole potential (0), C/(k) = (27r)^'^/^ J exp[ik-r] U{r) d^r, 
can be shown to be [pljT^], for the laser propagation and polarization as in Fig. |l|, 

^ ' ~ 2(2^)3/2e2c 3 ^ fc2 + _ q)2 ^ ^1 - q^ ^ kl + {ky + qf + kl ~ q^ ) ' ^' 

We proceed by adopting a cylindrically symmetric (about the axial z direction) variational ansatz for the density. 
For this to hold true requires the tight radial trapping to prevail over the (anisotropic) dipole-dipole forces in the 
radial direction. It is then reasonable to approximate the radial density by a Gaussian profile whose width is the 
variational parameter Wr'- n{r) = N (ttw^)"^ n^{z) exp [— (x^ -|- y'^)/wl.] , where N is the total number of atoms and 
n^{z) is the as yet unspecified axial density. Then, in momentum space, the density becomes 



m = (^"'(fc-) y-^iK + K)) ' (10) 

n^{kz) being the Fourier transform of the axial density n^{z). Assuming tight radial confinement with respect to 
the laser wavelength, Wr = Wr/Xh ^ 1, the radial integration in H^d can be approximately evaluated, so that the 
dipole-dipole energy reduces to a one dimensional functional along the axial direction 

Hdd = ^ j n\z)n'-{z')U'{z-z')dzdz' ^^-^ j n^-{k,)n^{-k,)U^-{k,) dk, , (11) 

where h{k) is the Fourier transform of the number density. In this expression the axial momentum has been scaled 
by the laser wavenumber kz = kz/q, and the one-dimensional (ID) axial potential is given by 

with 



4' ^ 2^V 2 / V 2 

1(1 - fc,^)(2.W.)^exp K^l~m^^r?\ , ^ (13) 
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where Ei is the exponential integral . 

We plot ?P(fc^) in Fig. | It tends to a positive constant value for large kz, as verified by the appropriate limit of 
(H) and (|l|): 

This means that the axial dipole-dipole potential in coordinate space will contain a repulsive delta-function contri- 
bution^ which can stabilize the trapped condensate in addition to (or instead of) the short-range s-wave scattering. 

Despite this effective short-range repulsion, the overall long-range behavior is dictated by the fcz — > limit of U^{kz), 
which is attractive. When the kz = Fourier component of the total ID-reduced potential (the sum of the s-wave 
scattering pseudo-potential and the laser-induced dipole-dipole interaction) becomes negative, the condensate size 
decreases drastically until it becomes self-trapped within the laser-wavelength. 
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FIG. 2. The dependence of the one dimensional axial potential U^{kz). The radial width has been set at Wr — l/(47r). 

From Fig. ^, the monotonic form of the axial potential suggests that the simplest ansatz for the axial variational 
density is also a Gaussian 



n^'ikz 



1 



exp 



(27rW^)2p 



(15) 



Even with this choice, however, the remaining integral in the calculation of H^d is difficult and will here be evaluated 
numerically. In the large N limit, the so-called Thomas-Fermi regime, i/ko may be neglected in comparison with the 
interaction energies and thus the energy i/xF = Hg + H^d + -ffho can be written in the dimensionless form 



Htf — 



1 



2'KWz{2TTWr)'^ 



{kz)Q{wr,kz)dkz 



1 



2Naq 



{2TrWr)^ (27rW^)^ 



(16) 



where ZJ^^ — \Jhj (mtOr) is the harmonic oscillator length due to a radial trap of frequency ujr, and similarly for 
^ho- Given these trapping frequencies, the real control parameter that allows one to play with the ratio of these 
dipole-dipole forces to the s-wave scattering is the dimensionless 'intensity' 



I = 



1 a^m 



(17) 



A minimization of H^f brings out the strong dependence of the condensate radius Wz and the weaker dependence of 
the Wr radius upon X as shown in Figures [s] and0. For X ^ 1.1 there is a drastic decrease of Wz and self-trapping in 
accordance with the discussion following Eq. (|14|). For X — 3/2 the system collapses and this can be attributed to 
the instability caused by the static r~^ part of the dipole-dipole potential (Q), as in ||l0|] . 
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FIG. 3. Condensate size as a function of the scaled laser intensity T. Solid line: w^, dashed line: Wr- The trap frequencies 
have been chosen so that at T = (i.e. without dipole-dipole forces) the condensate has an aspect ratio of 30:1, with a radial 
dimension oiwr ~ l/(47r). 




FIG. 4. Condensate size as a function of the scaled laser intensity J, as in fig. 3. here the trap frequencies are chosen 
such that at X = (i.e. without dipole-dipole forces) the condensate has an aspect ratio of 100:1, with a radial dimension of 
Wr = l/(47r). 

To conclude, we have predicted here a new feature of laser-irradiated cigar-shaped condensates: their strong com- 
pression along the soft axial direction. In an experiment this compression would provide a clear signature of laser- 
induced dipole-dipole forces at work. 
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